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The nature of dark matter is increasingly constrained by cosmological data. In this paper, we 
examine the implications of the Cosmic Microwave Background anisotropy limits on the density 
of cold dark matter under different theoretical assumptions and combinations of datasets. We 
infer the constraint S7 c[ j m /i 2 = 0.12 ± 0.04 (at 95% c.l.). The CDM models are compared with 
the shape of the linear matter power spectrum inferred from the 2dF galaxy redshift survey and 
with the rms mass fluctuations from recent local cluster observations. We found that a value 
of erg ~ 1 as suggested by recent cosmic shear data is not favoured by the CMB data alone 
nor by combined CMB+SN-Ia, CMB+HST or CMB+2dFGRS analyses. We also extrapolate 
our bounds on the rms linear mass fluctuations to sub-galactic scales and compare them with 
recent lensing constraints, finding agreement with the standard ACDM model. 

ations in spheres of 8/i _1 Mpc, erg- 

We extend and complement recent studies (see e.g. 
fl§], f||, H|) by studying the stability of the con- 
straints towards different theoretical assumptions such 
as, for example, a dark energy equation of state pa- 
rameter wq > —1 (quintessence) or an extra back- 
ground of relativistic particles (AN^M > 0). Fur- 
thermore, under the assumption of CDM, we extrap- 
olate the cosmological constraints into a measurement 
of the linear rms density fluctuations on sub-galactic 
mass scales. 



I. INTRODUCTION 

The new results on Cosmic Microwave Background 
(CMB) anisotropy from TOCO ( @, BOOMERanG ( 
II, H), MAXIMA ( [|), and DASI § experiments 
represent an extraordinary success for the standard 
cosmological model of structure formation based on 
Cold Dark Matter (CDM) and adiabatic primordial 
perturbations (see e.g. Jll|). Furthermore, early data 
releases from the 2dFGRS and SDSS galaxy redshift 
surveys ( JlTj, psf ) are living up to expectations and 
combined analysis of all these datasets are placing 
strong constraints on most cosmological parameters 
( & ©)■ 

However, even if theory and observations are in 
spectacular agreement on large ( }t IMpc) scales, var- 
ious discrepancies seem to be present on smaller (sub- 
galactic) scales. In particular, numerical simulations 
of CDM models yield an excess of small-scale power, 
producing, for example, more satellite galaxies than 
observed ( |20|) and cuspy galactic halos with exces- 
sive dark mass concentrations ( ||, Q), while on the 
smallest scales probed, gravitational lensing ( (!]], Q) 
prefers CDM substructure. Many authors have ad- 
dressed these issues and many solutions have been 
proposed, ranging from modifications of the proper- 
ties of the CDM particles to the introduction of new 
physics and to a refinement of the astrophysical pro- 
cesses involved in the simulations. 

One way to address these issues is by detection of 
dark matter. The key quantity of relevance to dark 
matter searches is the CDM density, which specifies 
the dark matter annihilation cross-section. In this pa- 
per, we use the latest CMB and deep redshift surveys 
data in order to simultaneously constrain the most 
important cosmological observables related to CDM: 
the CDM density £l c d m h 2 , the spectral tilt n$ of the 
primordial power spectrum and the rms mass fluctu- 



II. COSMOLOGICAL CONSTRAINTS ON 
CDM 

A. Method 



We compare the recent CMB observations with 
a set of models with 6 parameters sampled as fol- 
lows: Qcdmh 2 = uJcdm = 0.01, ...0.40, in steps of 
0.01; n b h 2 ee uu b = 0.001,. ..,0.040, in steps of 0.001; 
51a = 0.0, ...,1.0, in steps of 0.05 and such that 
Vt m = 0.1, ...,1.0, in steps of 0.05. The value of the 
Hubble constant is not an independent parameter, 
since: 



h = 



Ucdm + <^b 

i-n k -n A ' 



(i) 



We vary the spectral index of the primordial density 
perturbations within the range n s — 0.60, 1.40 (in 
steps of 0.02) and we rescale the amplitude of the fluc- 
tuations by a pre- factor Cio, in units of Ci OBE , with 
0.50 < Cio < 1.40. In order to test the stability of our 
results under variation of some of the theoretical as- 
sumptions, we also consider variations in the dark en- 
ergy equation of state parameter wq = —1.0, —0.5 
as expected in quintessence models (see e.g. p5j, 0|, 
fll3f ) and of an extra background of relativistic 
particles with AN„ — 0.3, 9 in steps of 0.3 (see [ pT~[ 
and references therein). 
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FIG. 1. Top panel: CMB anisotropics and CDM. Bot- 
tom panel: Allowed region for the matter power spectrum 
from CMB and from other cosmological observables ob- 
tained under the assumption of adiabatic CDM primor- 
dial fluctuations. The data from the 2dF redshift survey 
(Tegmark and Hamilton, 2002) is also plotted in the figure. 



The theoretical models are computed using a 
modified version of the publicly available CMB- 
fast code |3(|] and are compared with the recent 
BOOMERanG-98, DASI and MAXIMA-1 results. 
The power spectra from these experiments were es- 
timated in 19, 9 and 13 bins respectively, spanning 
the range 25 < I < 1150. For the DASI and 
MAXIMA-I experiments we use the publicly available 
correlation matrices and window functions. For the 
BOOMERanG experiment we assign a flat interpola- 
tion for the spectrum in each bin £(£+ l)Ci/2n — Cb, 
we approximate the signal Cb inside the bin to be 
a Gaussian variable and we consider ~ 10% cor- 
relations between the various bins. The likelihood 
for a given cosmological model is then defined by 
-21n£ = (C** - C b x )Mbb\C^, - Cjp) where M BB ' 
is the Gaussian curvature of the likelihood matrix at 
the peak. We consider 10%, 4% and 5% Gaussian dis- 
tributed calibration errors for the BOOMERanG-98 
[?], DASI and MAXIMA-1 [?] experiments respec- 
tively and we included the beam uncertainties by the 
analytical marginalization method presented in f22|j . 
We also include the COBE data using Lloyd Knox's 
RADPack packages. 
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FIG. 2. Constraints in the (fi m )°' 6 <78 — Q, c d m h 2 plane. 
The results of the 3 combined analysis CMB+HST, 
CMB+SN-Ia and CMB+2dFGRS are shown together with 
the 68% c.l. cluster constraints 



In addition to the CMB data we incorporate the 
real-space power spectrum of galaxies in the 2dF 100k 
galaxy redshift survey using the data and window 
functions of the analysis of Tegmark et al. ( |37j ) . 

To compute C 2dF , we evaluate pi = P(h), where 
P{k) is the theoretical matter power spectrum and 
ki are the 49 k- values of the measurements in [ B7j . 
Therefore we have -2ln£ 2dF = £\[P - {Wp) l ] 2 /dP 2 , 
where Pj and dPi are the measurements and corre- 
sponding error bars and W is the reported 27 x 49 
window matrix. We restrict the analysis to a range 
of scales where the fluctuations are assumed to be in 
the linear regime (k < 0.2/i _1 Mpc). When combin- 
ing with the CMB data, we marginalize over a bias b 
considered to be an additional free parameter. 

We attribute a likelihood to each value of uj c dm, 
ns, as by marginalizing over the nuisance parameters. 
We then define our 68% (95%), confidence levels to be 
where the integral of the likelihood is 0.16 (0.025) and 
0.84 (0.975) of the total value. 



CMB results and test for theoretical 
assumptions. 



The results of our analysis are shown in Table 1. 
In the first row, we restrict our analysis to CMB data 
with a combination of "weak priors": h — 0.65 ± 0.2, 
age t > lOGyr and tt m > 0.1. Since the value of 
as is degenerate with f2 m , in the table we report cr| 
defined as the value of erg with fl m = 0.30 ± 0.05. 
We see that, under the class of models considered, 
the CMB data suggest a ~ 2a detection of CDM, a 
scalar spectral index ns ~ 0.92 and a value of as ~ 
0.66 at fi m ~ 0.3. The reason why the present data 
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FIG. 3. Constraints in the Q m — erg plane. The results 
of the 3 combined analysis CMB+HST, CMB+SN-Ia and 
CMB+2dFGRS are shown together with the 68% c.l. con- 
straints from Viana et al. 2001 and Pierpaoli et al 2001. 



TABLE I. Parameter estimates from the CMB data sets 
and complementary data. Below the line, we restrict the 
parameter space to Q = 1. erg is the value of erg computed 
in the region Q m — 0.30 ± 0.05. Central values and 2a 
limits are found from the 50%, 2.5% and 97.5% integrals 
of the marginalized likelihood. 
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favour CDM is evident in Fig.l, top panel, where we 
plot the recent CMB data with the best fit purely 
baryonic (BDM) and CDM models. BDM models fail 
to reproduce the observed power at I > 700. The 
new CMB data provide independent support for the 
presence of non-baryonic dark matter in the universe. 

Before including information from complementary 
cosmological datasets, it is important to test the sta- 
bility of our CMB results by removing some of the 
theoretical assumptions used in the analysis. We re- 
strict our analysis to ACDM models. However an- 
other candidate that could possibly explain the ob- 
servations of an accelerating universe is a dynamical 
scalar "quintessence" field. The common character- 
istic of quintessence models is that their equation of 
state, wq = p/ p, varies with time and can be greater 
than wq = — 1, the value corresponding to the cosmo- 
logical constant. Adopting 'quintessence' instead of a 
cosmological constant does not change the results of 
our analysis, but increases the error bars. 

Another possibility is to consider an extra back- 
ground of relativistic particles (see e.g. 0), 
parametrized by a larger number of effective mass- 
less neutrinos Nyf. Since increasing N^* changes 
the epoch of equality, which is well determined by 
CMB observations ( UW), larger values of £l c d m h 2 are 
needed in order to compensate for this variation. The 
constraints on <r| are not greatly affected since the 
matter power spectrum is mainly sensitive to changes 
in the redshift of equality, which is kept constant via 
the CMB data. 

In the entire analysis, we assume negligible reioniza- 
tion and an optical depth r c ~ 0. This is in agreement 
with recent estimates of the redshift of reionization 
z re ~ 6 ± 1 (see e.g. IpH ). We have also removed this 
assumption. Due to the well known degeneracy with 
the scalar spectral index ns, the effect of including 



variations in r c leaves the CMB data in better agree- 
ment with higher values of ns and of erg. 

A background of gravity waves and/or of isocurva- 
ture modes can modify the theoretical CMB spectrum 
on large angular scales. In order to test for these hy- 
potheses, we repeat the analysis without the COBE 
data. This has the effect of relaxing our constraints. 

The high-£ part of the observed spectrum can be 
contaminated by different systematics (see e.g. p2|): 
beam reconstruction, detector noise, foregrounds. All 
3 different experiments can be affected by different 
systematics and the fact that the 3 datasets are in 
reasonable agreement suggests that the systematics 
are under control. However, we repeated the analysis 
removing the data points at £ > 650. As we can see, 
again, apart from an increase in the error bars, there 
is no significant difference. 

C. Comparison with complementary datasets. 

Since the CMB results are stable under variations of 
different theoretical assumptions, we can now assume 
that the A-CDM models are valid and investigate the 
effects of applying various prior probabilities and of 
incorporating complementary cosmological datasets. 
Including the gaussian prior 0.8f2 m — 0.60a = —0.2 ± 
0.1 from type la supernova luminosity distances 
or h = 0.71 ±0.07 from measurements with the Hubble 
Space Telescope |23) does not change our conclusions 
and improves the constraints on the 3 parameters. 

In Figure 1, bottom panel, we check for consistency 
of the 2dF data with the set of CDM models used 
in the CMB analysis by plotting a convolution of all 
the matter power spectra from the theoretical models 
in agreement with the CMB data, together with the 
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recent 2dF analysis of |57|. As one can see, the re- 
gion consistent with CMB alone is quite broad (due 
to the weak CMB constraint on f^) and contains the 
shape of the 2dF spectrum. Including other cosmolog- 
ical constraints from SN-Ia and HST shrinks the CMB 
constraint into a region consistent with the shape in- 
ferred from 2dF. This 'consistency' is reflected in the 
results of Table 1 : including information on the shape 
of the 2dF matter power spectrum improves our con- 
straints in a similar direction to the SN-Ia and HST 
priors. 

On similar scales, recent analyses of the local clus- 
ter number counts can be summarised as giving dif- 
ferent results for erg mainly due to systematics in the 
calibration between cluster virial mass and tempera- 



ture: a high value ~ Q^f&a 
ment with the results of 



and |3_ 

in the fl^fos ~ Qcdmh 2 plane 



0.55 ± 0.05 in agree- 
and a lower one, 
= 0.40 ± 0.05 following the analyses of §| 
It is interesting to plot these constraints 
We do this in Fig- 
ure 2, where we plot the 95% confidence level con- 
tour of the combined CMB+HST, CMB+SN-Ia and 
CMB+2dFGRS analyses (obtained again with the as- 
sumption of CDM) together with the high and low 
constraints on f^ 6 er 8 at 68% c.L. 

As we can see, a correlation appears from the 
present CMB data between these 2 quantities: 
namely, increasing fl c dmh 2 enhances the rms mass 
fluctuations. The independent constraint on fi^ erg 
from clusters can be used to break the degeneracy. 
Using the high constraint, we obtain a higher value 
for the density in CDM (see Table 1) with 0.17 > 
Qcmdh 2 > 0.11 at 95% c.L. Using the low value, 
the constraint becomes 0.15 > Vl cm dh 2 > 0.09, again 
at 95% c.L. The 2 results are consistent and favour 
CDM; however, when additional information such as 
SN-Ia and 2dF are included, the CMB tends to prefer 
the lower value. 

We further analyze this possible discrepancy in Fig- 
ure 3 where we plot constraints in the Q m — erg plane 
together with the 68% c.L results of j|| and of @. 
For values of Q m < 0.35, the high er 8 constraint is 
in slight disagreement with the CMB+2dF result at 
more than ler. 

On smaller (sub-galactic) scales, recent and strong 
constraints on the linear rms mass fluctuations cr(M) 
have been obtained through lensing measurements 
( 0). In Figure 3, we plot the predictions of the CDM 
models that are within the 95% CMB constraints in 
the er(_R) — R plane together with a region indicative 
of the constraint obtained from lensing. A more care- 
ful comparison and a better study of all the assump- 
tions and possible systematics should be done. How- 
ever, here we wish to point out that, even in our sim- 
ple analysis, the two contours overlap and the CDM 
paradigm seems consistent over a range of scales from 
lO^pch- 1 to lO^Mpctr 1 . 




R [Mpc h ] 

FIG. 4. The predictions of the CDM models that are 
within the 95% CMB constraints in the <j(R) — R plane 
together with a region indicative of the constraint obtained 
from lensing. 



III. DISCUSSION 

In this paper we have provided strong constraints 
on the amount of non-baryonic dark matter in the 
universe by comparing adiabatic inflationary mod- 
els with a set of cosmological observations. Com- 
bining CMB anisotropy measurements, high redshift 
supernovae observations, constraints on the Hubble 
parameter from the HST key project, and the mat- 
ter power spectrum data from the 2dFGRS, we find 
ttcdmh 2 = 0.11+g .eg, n s = 0.93±0.08, erg = 0.66±0.10 
at 95% confidence level. These results provide strong 
evidence for the presence of non-baryonic dark matter 
in a way independent of the local cluster abundance 
data. 

If the dark matter is the neutralino, a knowledge of 
Qcdmh 2 can set strong bounds on the parameter space 
of the simplest and most direct implementations of 
supersymmetry (see e.g. Q, [II]! , |ft3| ). 

The constraint on ^l c dnJ^ obtained here can be con- 
sidered stable under the removal of some of the the- 
oretical assumptions made in the analysis. Including 
quintessence, reionization and a large-scale CMB com- 
ponent such as expected from gravity waves (see e.g. 

33], J34j) or isocurvature CDM perturbations (see e.g. 

35]) can relax the constraints but does not change 
the conclusions. Models with £l c d m h 2 as large as 0.2 
are disfavoured in all cases except under the exotic 
hypothesis of an extra background of relativistic par- 
ticles. 

Depending on the external data we incorporate, the 
CMB constraints in the os&lf — ^cdmh 2 plane can be 
1 — 2cr lower than most of the determinations inferred 
from the local cluster X-ray temperature function (se e 
e^ [fpL j25|], [||]) and cosmic shear data ( |f], 
eJL IH5|), while in better agreement with the 'new' 
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analyses of p9] and Q. However, we showed that 
use of the high or low priors on as has only a marginal 
effect on the constraints on f2 ct ; m /i 2 . 

With reference to previous analyses, we are in agree- 
ment with the recent paper by Lahav et al. [Q, al- 
though our value of <7g appears to be slightly lower. 
However we allow for variations in the spectral index 
ns while in |?8| this parameter has been set tons = f . 
We also included the COBE data using the lognormal 
approximation as in jIo| and the 2dF data is taken 
from the independent analysis of Q . 

The low values of as and ns can possibly alleviate 
some of the problems of CDM on subgalactic scales. 
We will investigate this in a forthcoming paper pH| . 
However, we also find that the models compatible with 
the CMB can satisfy the recent lensing constraints. 
This result suggests that a solution of the subgalactic 
CDM problem can most likely be obtained by a re- 
finement of the astrophysical processes involved in the 
numerical simulations, rather than by an ad hoc mod- 
ification of the dark matter properties. However, if fu- 
ture cluster temperature or cosmic shear analyses were 
to converge towards a higher as value, then this could 
lead to a possible discrepancy with the CMB+2dF re- 
sult. It will be the task of future experiments and 
analysis to verify this interesting result. 
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